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We observe asymmetric transition rates between Zeeman levels (spin-flips) of magnetically trapped 
atoms. The asymmetry strongly depends on the spectral shape of an applied noise. This effect 
follows from the interplay between the internal states of the atoms and their external degrees of 
freedom, where different trapped levels experience different potentials. Such insight may prove useful 
for controlling atomic states by the introduction of noise, as well as provide a better understanding 
of the effect of noise on the coherent operation of quantum systems. 
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PACS numbers: 37.10.Gh, 32.70.Cs, 05.40.-a, 67.85.-d 



Coupling between internal and external (motional) de- 
grees of freedom plays a major role in cooling atoms, 
ions and molecules [1- 3] and in manipulating their quan- 
tum states, e.g., in logic gates for quantum computation 
0-0|- The state of these particles is usually controlled 
by monochromatic (or transform limited) electromag- 
netic fields, while incoherent fluctuations (noise) must 
be suppressed in order to prevent decoherence, heating 
and loss The latter are usually studied under the 

"white noise" assumption. Little is known about what 
happens between the monochromatic and white noise 
limits ("colored noise") with respect to control and hin- 
dering effects. 

To understand how the noise spectrum can affect the 
rate of transitions between internal states, consider a 
system of two levels representing electronic configura- 
tions of an atom, ion or molecule, coupled to external 
(translational, rotational or vibrational) degrees of free- 
dom and to a weak homogeneous field inducing tran- 
sitions between the internal levels. If this field im- 
poses a monochromatic perturbation \e~ lut , then the 
rate of transitions from an initial state \i) to a fi- 
nal state |/) is given by Fermi's golden rule 7i_>./ = 
(2ir/S) Efe KW, k)\ 2 S(Ei - Ef, k + M, wh ere k repre- 
sents the quantum numbers of the external degrees of 
freedom of |/). The transition rates are proportional to 
the density of states in their respective final level, and are 
therefore asymmetric between two levels that experience 
different external potentials. At the other extreme, one 
has a fluctuating random field Xg(t), where g(t) has a 
spectral density S(u) = J dt'e iu ^- t "> {g(t)g(t')) which is 
flat over a large bandwidth ( "white noise" ) . This now al- 
lows transitions to all k states whose completeness yields 
7i-»-/ = 5|(i|A|/)| 2 such that the external degrees of free- 
dom decouple from the transition dynamics and the tran- 
sition rates between the two levels become symmetric. 

Numerous quantum systems proposed for applications 
such as quantum information processing are based on 
particles trapped in external potentials. A major limita- 
tion of these systems is uncontrolled noise-induced transi- 
tions between the internal states used for the application 



(e.g. qubits) or from the latter to other internal states. 
Typically, different states may experience different po- 
tentials, and at times the potential difference is made 
significant as part of the application itself (e.g. quantum 
logic gates). As conventional environments usually con- 
tain background or technical noise which is not white, 
it is important to understand the interplay between the 
internal and external degrees of freedom under these con- 
ditions. 

In this work we study the dynamics of transitions be- 
tween two internal atomic states with different external 
potentials in the form of two magnetically trapped Zee- 
man levels (\F, ttif) — |2, 2) and |2, 1)) in the presence of 
colored noise which is neither monochromatic nor white. 
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FIG. 1. (a) The measured ratio R(t) between the number of 
atoms in iuf = 1 and the number of atoms in both trapped 
levels. Different curves correspond to different values of A/, 
the detuning between the noise peak frequency fo and the 
Zeeman splitting Ei 2 /h between the two trapped levels at the 
magnetic field minimum. The measured R(t) goes beyond the 
band < R < 1/2 representing asymptotic values of R(t) in 
a model where the external degrees of freedom are decoupled 
from the transition dynamics (the dashed line at R = 1/3 
is for white noise). Curves are fits to the empirical form of 
Eq. (fT|. Data variance gives an error of ±0.01 (not shown), 
(b) The noise spectrum we introduce into the system relative 
to fo — 18 MHz (red line is the fit used in our calculations). 
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We experimentally observe that the relative transition 
rates between the levels strongly depend on the spectral 
shape of the noise. While a clear understanding of this 
dependence is important in order to find effective ways 
to combat uncontrolled noise, we demonstrate that this 
dependence also allows steering the transitions in the de- 
sired direction by utilizing engineered noise. 

We start our experiment with ~ 7 x fO 4 87 Rb atoms 
evaporatively cooled to ~ 350 nK. This leaves most of 
the atoms in the |2, 2) Zeeman level with a residual frac- 
tion of ~9% in |2, 1). The magnetic trap is produced 
250 /im from the surface of an atom chip by a right- 
angle Z-shaped wire ll| and a homogeneous magnetic 
field generated by distant coils. This creates a harmonic 
trap with an axial frequency ui x — 2tt x 10 Hz (calculated) 
and radial frequencies ui y = lu z = 2n x 96 Hz (measured) 
for the |2,2) level (a factor of 1/V2 less for the |2, 1) 
level). In order to study the effect of the noise spectrum 
we introduce controlled magnetic noise with the same 
function generator (Agilent 33250A) and antenna used 
for the evaporative cooling. To suppress uncontrolled 
noise, we increase the Zeeman splitting to 18MHz. At 
this range of spin-flip transition frequencies the uncon- 
trolled noise spectrum is quite small and flat. The pop- 
ulations Ni and N 2 of the trapped levels nip — 1 and 
nip = 2, respectively, are determined by a Stern- Gerlach 
measurement [liq . 

Fig. Q}l presents the ratio R(t) = Ni/(Nj. + N 2 ) as a 
function of the time for which the noise is applied, for 
different detuning A/ between the noise peak frequency 
fo and the Zeeman splitting E^/h of the two trapped 
levels at the magnetic field minimum. Fig. [TJd presents 
the spectrum of the noise we introduce. One may observe 
a large difference in the evolution of R(t) between red- 
and blue-detuned noise. When the noise is red-detuned 
(A/ < 0), the nip = 1 level is populated while when 
the noise is blue-detuned (A/ > 0) the mp = 1 level is 
depleted. The solid curves are fits to the empirical form 



R(t) = Roo + {Rq - Roo)e 
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(1) 



which represents an exponential convergence from the 
initial value i?o = R(t = 0) to an asymptotic value 
Roo = R(t — > oo). These results show that the rela- 
tive transition rates between the levels strongly depend 
on the detuning of the noise and significantly differ from 
the expected rates in the case of white noise (dashed line 
in Fig. [Tk) . The value of 7 in Eq. (TT|) strongly depends 
on the intensity of the noise at the relevant frequency. 
At resonance (A/ = 0) it is roughly hundreds of Hz, 
depending on the exact noise amplitude used. 

The transitions between the atomic levels are caused 
by the coupling A = —gpiisjP • B no i SO between the mag- 
netic field noise B no isc and the magnetic moment of the 
atom. Here /ib is the Bohr magneton, F is the an- 
gular momentum operator and gF is the Lande factor. 
This coupling allows only transitions with spin change 
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FIG. 2. (a) The dependence of [Eq. Q] on /? for var- 
ious values of a. (b) Asymptotic values Roa taken from ex- 
perimental curves as in Fig. [T] as a function of A/, com- 
pared with theory [Eqs. ([11)-©] with no fitting parameters 
(T — 0.5 — 1.5 fJ, K). Error bars are rms of data variance and 
mean fit error in Fig. [T] except for points marked in red, for 
which only one measurement is available and the error is the 
average of the whole data set. The horizontal error is the 
uncertainty in the magnetic field minimum. 



jAm^l = f, which implies the following rate equations 
for the populations of the two trapped levels 



Ni = -(71^2 + 7i->o)-Wi + 72^1^2 

N 2 = 71^2^1 ~ 72^1^2- 



(2) 



We estimate the escape time of an atom in the untrapped 
level mp = from the trapping region, due to thermal 
velocity and gravitational acceleration, to be ~ 10 ms. 
For I A/| > 20kHz (±10 kHz due to the variance in 7) 
the atom is lost from the trap before it can make a tran- 
sition back to mj? = 1 and the transition 7o->i may be 
safely neglected. In the limit of white noise, where the 
motional and internal degrees of freedom are decoupled, 
we expect the ratio between the transition rates to be de- 
termined by the matrix elements of the angular momen- 
tum operator, such that a = 71-^0/72^1 = 3/2, while 
P = 71^2/72^ 
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13(. In the following we show 
that colored noise dictates different values for a and j3. 

Next we connect the ratio R(t) to a and (3 and focus 
on the asymptotic value Roo . Under the assumption that 
the transition rates are time-independent, Eqs. yield 
the solution 



R(t) 



Nt(t) 



Rn 



N 1 (t)+N 2 (t) 1 + aRooCe-^ 



(3) 



where C = (R - Roo) /(I - ai? oo i? ) and 7 = (l/Roo - 
ai?oo)72^i- Note that Eq. © reduces to Eq. ([T]) when 
a = 0, and even when a ^> 1 Eq. ([1]) is a fair approxima- 
tion. The asymptotic value in Eq. ([3]) is given by 



R? 



l+a + 0- ^-4a + (1 + a + f3) 2 
2a 



(4) 
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This value depends solely on the ratios a (determining 
the relative depletion rate to the roj? = level) and (3 
(determining the asymmetry in the transition rate be- 
tween the two trapped levels) , and is independent of the 
overall noise intensity or the initial condition, Rq. In the 
white noise limit, where a — 3/2 and /3 = 1, Roo = 1/3 
(dashed line in Fig.[T^) and Roo < 1/2 for any value of a 
(/3 = 1 band in Fig. [Tk). The dependence of Eq. (j4|) on 
j3 for various values of a is presented in Fig. [2ji. 

In Fig. [2Jd we plot the measured value of Roo for dif- 
ferent values of the center frequency of the applied noise. 
The band and the three curves represent the theoretical 
prediction for Roo, as described below. 

In order to gain a qualitative understanding of the re- 
sults, we first present a simple semiclassical ID model 
of two thermal distributions of atoms trapped in po- 
tentials Vj(x) — hmjMuJiX 2 , where M is the atomic 
mass, u>i is the trapping frequency of the mp = 1 level 
and 77i j = mp is either 1 or 2. Each atomic distribu- 
tion is represented in Fig. [3] by a typical atom at an 



average position dj 



y/k B T/mjMwl. As 



the wavelength of the applied noise is much larger than 
the typical size of the system, we may assume that mo- 
mentum is conserved locally during the transition (the 
recoil is negligible). The transition of a typical atom 
from mp = 2 to mp = 1 requires a photon energy 
E 2 ^i = E° 12 + V 2 (d 2 ) - Vi(d 2 ) = E° 12 + \k B T and re- 
duces the energy of the atom relative to the trap bottom 
by \k B T . A transition of a typical atom in mp = 1 
to m p — 2 at d\ = \/2d 2 requires a photon energy 
Ei^ 2 = E\ 2 + V 2 (di) - Vi(di) = E\ 2 + \k B T and in- 
creases the energy of the atom by ^k B T. 

The relative transition rates depend on the number of 
photons with the two energies E\^2 > ^2->i ■ If the 
noise intensity increases with frequency in a wide enough 
band (~ k B T) around / = E\ 2 jh (see blue arrows in 
Fig. A/ > 0), the transition 1 —> 2 is preferred and 
the population of mp = 1 is depleted, as in the plateau 
on the right hand side of Fig. [5J3. If, by contrast, the in- 
tensity decreases in this band (see red arrows in Fig. |3p, 
A/ < 0), the transition 2 — > 1 is preferred and the popu- 
lation of mp — 1 becomes dominant, as in the plateau on 
the left-hand side. Due to these effects, even when the 
noise intensity decreases by a few orders of magnitude 
(side peaks of the noise), the atoms are still sensitive to 
the details of the spectrum, as can be seen in Fig. [5p. 
This may prove to be a useful tool for characterizing noise 
features. Note also that the temperature determines the 
width of the spectral region the cloud samples, hence 
colder clouds are more sensitive to the fine details of the 
noise (e.g. 1 /jK gives a resolution of 20 kHz). This effect 
is most noticeable at A/ ss 0.7 MHz (Fig. O3) where the 
temperature band becomes much wider. 

To quantitatively explain the results, the ratios a and 
/3 are calculated using a semiclassical expression for the 




position d,2 d\ 



FIG. 3. (a) A simplified model showing the asymmetry in 
the transitions between the two trapped levels. We assume 
that the atoms are in thermal equilibrium at temperature T, 
in both levels. We also assume momentum to be unchanged 
during the transition as the photon recoil is negligible, and 
this is why the length of the arrows corresponds to the differ- 
ence between the potentials. We plot the transitions at the 
mean atom distances from the trap center, dj = (x 2 )j. A 
spin flip from mp = 2 to uif = 1 requires photon energy 
-E2-H, which is smaller than the photon energy £1^2 needed 
for the reverse transition, (b) The two transitions sample 
different frequencies of the noise spectrum (plotted twice for 
opposite detuning), giving rise to the asymmetric transition 
rates. The blue and red pairs of arrows point to the sampled 
noise intensities for each detuning. 



transition rate which follows from Fermi's golden rule [Tj 



du)Tif(u>) 



G? 3 pd 3 i 



■Pi(r,p)J(AVi/(r)-fiu) 



where T if { U ) = Ei=»,, K*1^l/> »u 



(27T/7) 3 

(5) 
S 3 i be- 
ing the spectral density of the magnetic fluctuations [1, 
13] , and Pi(r, p) is a normalized phase space distribu- 
tion for level i. AV if = E? f + V/ - V i: where E$ f is 
the Zeeman splitting calculated using the Breit-Rabi for- 
mula, predicting a difference of h x 95 kHz between E® 2 
and E$ v Vj = |mjMS fclA! wL^ + Mgz, where u lk 
is the trapping frequency of the mp = 1 level along the 
A;th axis (frequencies given above), rtij (j — i, /) are the 
indices of the initial and final Zeeman levels, and g is 
the gravitational acceleration which causes a shift of the 
trap centers relative to each other. Assuming a Maxwell- 
Boltzmann distribution, Pi cx e~ p / 2mk BT e -Vi/k B T w ^ n 
a constant temperature T for both levels, Eq. (|5|) reduces 
to an integral over dimensionless space coordinates 



3/2 



7i->/ = 4- 



dgg 2 c -(WW/m,) sinh (2??g) 



2i]q 



xr lf (Lu = (E° f + q 2 k B T)/h), 



(6) 



where r\ = (g/uji z ) y/ M /2k B T . Trap frequency variations 
due to the non-linear Zeeman effect arc neglected. 

The theoretical curves are calculated utilizing a fit to 
the independently measured spectrum of the noise (red 
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curve in Fig. [T]d): a Lorentzian of 1kHz FWHM mul- 
tiplied by a Gaussian with a = 150 kHz for the center 
peak, two identical Gaussians for the two side peaks, and 
a constant for the background white noise. Finer details 
of the noise are ignored. The band in Fig. [5b represents 
the different theoretical values of Roo for cloud tempera- 
tures ranging from 0.5 /iK to 1.5 /iK. The 1 /iK wide band 
accounts for heating due to the transition process (\UbT 
per 2—^1—^2 cycle in Fig. 3) and due to uncontrolled 
background noise. 

To further verify our model and to demonstrate control 
by noise, we conducted the following experiments: first, 
we introduce red-detuned noise (A/ = —0.2 MHz). At 
this frequency the system has a steady state of Roo s» 
0.7, as can be seen in Fig. [5Jd. When most of the atoms 
are in the mp — 1 level, we rapidly change the center 
frequency to A/ = +0.4 MHz. In the new steady state 
all atoms should be in the mp — 2 level, exactly as we 
have observed (Fig. 0^,) . 

It is interesting to consider our model under the 
two limits noted previously, namely, white noise and 
monochromatic radiation. For white noise our theory 
predicts Roo = 1/3, which is confirmed by our measure- 
ment presented in Fig. The prediction of our model 
for the other limit may be estimated by introducing into 
Eq. ^ Fif(u}) oc 5(u — Wjioise)' It follows that transitions 
occur only for Af > and that the ratio (3 becomes 



0n 



= 2- 3 / 2 



exp 



2%hAf - Mg 2 /4uj 2 z 



(7) 



The factor 2~ 3 / 2 is the 3D ratio between the density of 
states in the two levels. The second factor in the ex- 
ponent is due to the different gravitational shift of the 
atomic potential minima. When the transition 1 — > is 
far detuned from the transition 2 — > 1 (a — > 0) due to 
the non-linear Zeeman shift, lim Q ,_>.o Roc = 1/(1 + /3) ~ 
0.92 (1 ijK), 0.85 (2 fjK), relative to a value of Roo « 0.73 
expected from the density of states ratio alone. This ex- 
plains the high theoretical values in Fig. [^b near reso- 
nance, where the noise is dominated by a sharp peak. We 
attribute the fact that Roo < 0.8 in all our near-resonance 
measurements to the fast heating rates which were ob- 
served independently. Note also that in this region our 
assumption regarding thermal equilibrium is invalid. At 
| A/| < 150 ± 50kHz, the collision rate (estimated to 
be ~ 0.5 Hz) is slower than the spin flip transition rates 
estimated by 7. Nevertheless, asymmetry exists in this 
regime as well, which may indicate that thermalization is 
not the fundamental source of the observed asymmetry. 
Indeed, in Eq. ([SJ we have used the Maxwell-Boltzmann 
distribution as a mere simplifying assumption and we 
believe that a more elaborate model would reach similar 
results without invoking thermal equilibrium. 

To summarize, we have observed that the spectral 
shape of noise determines the relative transition rates be- 
tween internal levels of trapped neutral atoms, due to the 




FIG. 4. (a) Demonstration of control by noise. We use red- 
detuned noise (A/ = —0.2 MHz) for 200 ms transferring most 
of the atoms into level rriF — 1 (as predicted in Fig. db). 
We then (t = 0) apply blue-detuned noise (A/ = +0.4 MHz) 
which transfers all the atoms into level ttif = 2, as predicted 
by theory (Fig. [2b). (b) Measured ratio R(t) for applied white 
noise. The measured asymptotic value of Roo = 0.344 + 0.004 
is very close to the predicted Roo = 1/3. Errors are estimated 
as in Fig. [TJ 



interplay of these levels with the external degrees of free- 
dom. As non-white noise spectra are common in back- 
ground and technical noise (e.g. in the radio-frequency 
regime), and may be very significant in some types of 
apparatuses such as atom chips [15j, these findings may 
serve to better understand how noise couples to quantum 
systems based on trapped atoms. Future studies should 
investigate the effects of colored noise on the coherence 
in systems of trapped particles with few levels, and the 
possibility that colored noise may be used for the sup- 
pression of decoherence, as proposed in Ref. 



16] 



RF thanks Yoseph (Joe) Imry for being a profound 
mentor. SM would like to thank Julien Chabe for helpful 
discussions. We thank A. Aharony, O. Entin-Wohlman, 
B. Horovitz and M. Schechter for stimulating discussions. 
We acknowledge support from the German-Israeli fund 
(GIF) and from the Israeli Science foundation (ISF). 



Am. 



[2] 

[3] 

[5] 
[6] 
[7] 

[8] 



machlufs@bgu.ac.il 

Present address: Center for Cold Matter, Imperial Col- 
lege London, UK 

J. Dalibard and C. Cohen- Tannoudji, J. Opt. Soc. 
B 6, 2023 (1989). 

Phys. Rev. Lett. 75, 4011 (1995). 
., Science 321, 1062 (2008). 
P. Zoller, |Phys. Rev. Lettl, 74, 



C. Monroe et al. 
J. G. Danzel et 1 
J. I. Cirac and 
(1995). 

T. Calarco et al, Phys. Rev. A 61, 022304 (2000). 

O. Mandel et al, Nature 425, 937 (2003). 

M. F. Riedel et al, Nature 464, 1170 (2010). 

C. Henkel, S. Potting, and M. Wilkens, Appl. Phys 

69, 379 (1999). 



4091 



■5 



[9] M. P. A. Jones et al, Phys. Rev. Lett. 91, 080401 (2003). 
[10] Y.-J. Lin et al, Phys. Rev. Lett. 92, 050404 (2004). 
[11] S. Wildermuthet al, Phys. Rev. A 69, 030901(R) (2004). 
[12] A. E. Leanhardtet al, Phys. Rev. Lett. 90, 100404 
(2003). 

[13] V. Dikovsky et al, Eur. Phys. J. D 35, 87 (2005). 



[14] Y. Japha and B. Segev, Phys. Rev. A 65, 063411 (2002). 
[15] R. Folman et al, Adv. At. Mol. Opt. Phys. 48, 263 
(2002). 

[16] G. Gordon, N. Erez and G. Kurizki, J. Phys. B 40, S75 
(2007); Y. Khodorkovsky, G. Kurizki and A. Vardi, Phys. 
Rev. Lett. 100, 220403 (2008); N. Bar-Gill et al, Phys. 
Rev. A 80, 053613 (2009). 



